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ABSTRACT 

We use observations from the Giant Metrewave Radio Telescope (GMRT) to mea- 
■ sure the atomic hydrogen gas content of star-forming galaxies at z = 0.24 (i.e. a look- 

j — | backtime of ^3 Gyr) . The sample of galaxies studied were selected from Ha-emitting 

, field galaxies detected in a narrow-band imaging survey with the Subaru Telescope. 

The Anglo- Australian Telescope was used to obtain precise optical redshifts for these 
galaxies. We then coadded the HI 21 cm emission signal for all the galaxies within the 
GMRT spectral line data cube. From the coadded signal of 121 galaxies, we measure 
an average atomic hydrogen gas mass of (2.26 ± 0.90) x 10 9 Mq. We translate this HI 
\ signal into a cosmic density of neutral gas at z = 0.24 of f2 gas = (0.91 ± 0.42) x 10~ 3 . 

This is the current highest redshift at which f2 gas has been constrained from 21 cm 
emission and our value is consistent with that estimated from damped Lya systems 
around this redshift. We also find that the correlations between the Ha luminosity 
and the radio continuum luminosity and between the star formation rate and the HI 
' gas content in star-forming galaxies at z = 0.24 are consistent with the correlations 

^ . found at z = 0. These two results suggest that the star formation mechanisms in field 

galaxies ~3 Gyr ago were not substantially different from the present, even though 
the star formation rate is 3 times higher. 

Key words: galaxies: evolution, galaxies: ISM, radio lines: galaxies, radio continuum: 
galaxies 



1 INTRODUCTION 

The cosmic star formation rate density is known to increase 
by an order of magnitude between the present time and z ~ 1 
(Lilly et al. 1996; Madau et al. 1996; Hopkins 2004). How- 
ever, the cosmic density of neutral gas is poorly observation- 
ally constrained over this redshift range. Consequently, we 
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are missing a key element in our understanding of galaxy 
evolution, that of the fuel supply for new stars. 

In the local universe (z ~ 0) the quantity of atomic 
hydrogen gas in galaxies has been measured with good pre- 
cision using HI 21 cm emission surveys (Briggs 1990; Zwaan 
et al. 1997, 2005). Observations of damped Lya absorption 
in QSO spectra have proven very effective at measuring 
the density of atomic hydrogen at high redshifts, z > 1.5 
(Storrie-Lombardi & Wolfe 2000; Prochaska et al. 2005). 
However, it is difficult to constrain f2 gas from observations of 



2 Lah et al. 



damped Lya absorbers at intermediate redshifts. The diffi- 
culties are: (i) that observations have to be done from space, 
since at these redshifts the Lya line is not accessible by 
ground-based telescopes, and (ii) that the average number 
of damped Lya absorbers per unit redshift is so low that 
an accurate estimate of fi gas requires impractical amounts 
of observing time. At lower redshifts (z < 1.5), Rao et al. 
(2006) have made measurements of the neutral gas density 
using Hubble Space Telescope (HST) observations of clouds 
that were preselected to have strong Mgll absorption. How- 
ever, this is a biased sample of objects compared to the ran- 
dom lines of sight probed by damped Lya systems at high 
redshift. In addition, the identification of the absorbing sys- 
tem in damped Lya measurements is problematic, as the 
bright QSO light makes imaging or spectroscopy difficult of 
the fainter galaxies hosting the absorber. 

Using radio observations of HI 21 cm emission to con- 
strain fig as at z > is difficult. The problem is that the flux 
from most individual galaxies at these redshifts falls below 
the detection limit of current radio telescopes for reasonable 
observing times. Our method for bringing the measurement 
of the HI 21 cm emission within reach is to coadd the signal 
from multiple galaxies using their observed optical positions 
and redshifts. This coadding technique for measuring atomic 
hydrogen gas has been used previously by Zwaan (2000) in 
galaxy cluster Abell 2218 at z = 0.18 and Chengalur et al. 
(2001) in galaxy cluster Abell 3128 at z = 0.06. We expand 
this technique to field galaxies with active star formation 
selected by their Ha emission. With this technique, we can 
directly relate gas content, star formation rate and stellar 
mass (from continuum luminosities) in the same sample of 
objects. 

The structure of this paper is as follows. Section 2 de- 
tails the optical and radio observations and data reduc- 
tion. Section 3 comments on the relationship between the 
1.4 GHz radio continuum and Ha luminosities at z = 0.24. 
Section 4 details our main results, the HI content of star- 
forming galaxies at z = 0.24. Finally, Section 5 presents a 
brief summary and discussion. 

We adopt the consensus cosmological parameters of 
Q,\ = 0.7, Q M = 0.3 and H = 70 kms -1 Mpc -1 throughout 
the paper. 



2 THE DATA 

2.1 The Optical Data 

The optical galaxies we are examining come from a deep 
narrow-band imaging survey for Ha emission galaxies at 
z = 0.24 (Fujita et al. 2003). This survey was performed 
using Suprime-Cam on the Subaru Telescope covering a 
field 34' x 27' in size. Deep observations were taken with 
broad passband filters B, Re, Ic, z' and a narrow pass- 
band filter NB816 centred at 8150 A (corresponding to 
the Ha emission line at z = 0.24). After applying vari- 
ous selection criteria to remove contaminants, Fujita et al. 
found a total of 348 Ha-emitting galaxies with equivalent 
widths greater than 12 A (hereafter these galaxies will be 
referred to as the 'Fujita galaxies'). From these galaxies, 
an extinction-corrected luminosity function for Ha emitting 
galaxies at z = 0.24 was derived. Using this luminosity func- 
tion, Fujita et al. computed a star formation rate density of 



0.036tomt M Q yr _1 Mpc" 3 . This is a factor of 3 higher than 
the local star formation rate density (Gallego et al. 1995; 
Condon et al. 2002; Gallego et al. 2002). Their measured 
star formation rate density is in good agreement with other 
measurements made at and around z = 0.24 (Hogg et al. 
1998; Tresse & Maddox 1998; Haarsma et al. 2000; Pascual 
et al. 2001; Georgakakis et al. 2003). This good agreement 
indicates that there is nothing particularly unusual about 
this field; it is not biased by containing an excess or defi- 
ciency of galaxies that have unusually high or unusually low 
star formation rates. 

We observed the 348 Fujita galaxies with the 2dF and 
AAOmega spectrographs at the Anglo Australian Telescope 
to obtain precise spectroscopic redshifts. Observations were 
carried out with 2dF on 6-10 March 2005, and an additional 
two hours of AAOmega service time obtained on 20 April 
2006. The data were reduced using the Anglo- Australian 
Observatory's 2dfdr data reduction software . As the tar- 
gets were selected for Ha emission, we confirmed their red- 
shift primarily using the emission lines of [OIIJA3727, H/3, 
[OIII]A4959 and [OIIIJA5007. We obtained secure redshifts 
for 154 galaxies with 0.2150 < z < 0.2641 and a maximum 
error in the redshift of 70 kms -1 . The break down of the 
redshifts by Ha luminosity is: 47 redshifts for the 50 galax- 
ies with L(Ha) > 40.9 log ergs -1 , 48 redshifts for the 58 
galaxies with 40.4 < L(Ha) ^ 40.9, and 42 redshifts for the 
240 galaxies with L(Ha) 40.4. 

2.2 The Radio Data 

Radio observations of this field were carried out on 26-28 
December 2003 and 7-9 & 24-26 January 2004 using the Gi- 
ant Metrewave Radio Telescope (GMRT) located near Pune 
in India. A total of 80.5 hours of telescope time was used 
with ~40 useful hours of on-source integration after the 
removal of slewing time, flux and phase calibrator scans, 
and the flagging of bad data. At most only 28 of the 30 
GMRT antennas were providing useful data on any one 
day. The total observing bandwidth of 32 MHz was split 
up into two 16 MHz- wide sidebands covering the frequency 
range from 1034 MHz to 1166 MHz (i.e. a redshift range 
0.218 < z < 0.253 for redshifted 21 cm HI emission). The 
pointing centre of the observations was R.A. 10 h 44 m 33.0 s 
Dec. -01°21'02" J2000 (see Figure 1). There were two po- 
larisations and 128 spectral channels per sideband, giving a 
channel spacing of 0.125 MHz (32.6 kms -1 ). Primary flux 
calibration was done using periodic observations of 3C286, 
which has a flux density at 1150 MHz of 16.24 Jy. Phase cal- 
ibration was done using scans on the VLA calibrator sources 
0943-083 and 1130-148, for which our observations give flux 
densities of 3.137±0.011 Jy and 5.163±0.044 Jy respectively. 

The data reduction was primarily done using AIPS and 
was partially automated using Perl scripts to interface with 
AIPS. Each sideband of data were processed separately. The 
initial assessment of data quality and the preliminary cali- 
bration was done using the visibilities in a single channel. 
Flux and phase calibration was determined using the band- 
pass calibration task BPASS, and the solution was inter- 
Anglo-Australian Observatory website 
http: / / www.aao.gov.au/AAO / 2df/aaomcga/aaomcga.html 
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Figure 1. This figure shows the positions of the optical galaxies 
as well as the pointing and beam size of the GMRT observations. 
The 348 Fujita galaxies are the small points. The large hexagonal 
point is the strong radio source J104420-011146 on which self 
calibration was performed. The small, bold, open, central circle 
is the GMRT pointing of the observations. The unbroken circle is 
the primary beam diameter of 29.2' and the broken circle is the 
10% level of the beam. The overlapping square grid shows the 16 
facets used to tile the sky in the radio imaging. The centre of this 
grid pattern was created with an offset from the GMRT pointing 
to ensure that all Fujita galaxies lie within the tiled region. 

polated to the data. No normalisation was done before de- 
termining the solutions. The regular phase calibrator scans 
throughout the observations allow us to correct for any time 
variability of the bandpass shape. Since, over the time span 
of the observations, there was no significant variation in the 
measured flux of the phase calibrators, the measured fluxes 
for each phase calibrator were combined into a single aver- 
age value. In the final data reduction, this average value was 
used as the calibrated flux for the phase calibrators. 

After a period of testing and refinement, we adopted 
an imaging pixel size of 0.7" and a robustness value of 0. 
For primary beam correction, we assumed a Gaussian beam 
with FWHM of 29.2' at 1150 MHz 2 . 

The brightest radio continuum source in the field is 
J104420-011146. It lies ~9.8' from the phase centre and its 
measured integrated flux (at 1150 MHz) is 303.66T0.05 mjy. 
The integrated flux listed in the FIRST radio catalog (at 
1.4 GHz) for this source is 276.83T0.15 mjy. J 104420-01 1146 
is an order of magnitude brighter than the next brightest ra- 
dio source in our field. Self-calibration of the data was done 
using this source. 

When making continuum images only channels 11-110 
(out of the 128 in each sideband) were used. In order to avoid 
bandwidth smearing, the visibilities were averaged into a 
new data set consisting of ten channels, each of which was 

2 National Centre for Radio Astrophysics website 
http: / / www.ncra.tifr.res.in/~ngk/primarybeam/beam 



Figure 2. This figure shows a grey scale image of one of the 
radio continuum facets (see the text for details). This facet is just 
south of the facet with the strong continuum source J 104420- 
011146 (see Figure 1). The residuals from this bright source run 
north-south and can be seen on the left hand side of the image. 
The RMS in the continuum images is 15 fijy, but along this strip 
of residual structure, the RMS increases to ~28 fijy (this is the 
worst residual remaining in the any of the continuum facets). 
The grey scale used in the image is a square root scale which 
varies from ~2.5 mjy for black and ~-120 /jJy for white. The 
brightest source in this image is in the top right hand corner and 
has extended structure. It been identified as the nearly face-on 
spiral galaxy, UGC 05849 at z = 0.026. The peak flux from the 
galaxy is 8.76 mjy/beam and the total flux density is 42.3 ± 
0.4 mjy, which is spread over a diameter of ~60". 



the average of 10 of the original channels. This new ten chan- 
nel u-v data file was turned into a single channel continuum 
image in the AIPS task IMAGR, combining the channels 
using a frequency- dependent primary beam correction based 
on the antenna size of 45 m. 

The self-calibration of GMRT data sometimes does not 
converge quickly, probably as a consequence of the GMRT 
hybrid configuration. Sources that lack coherence in the syn- 
thesis image due to phase errors tend to remain defocused. 
To fix this problem, slightly extended sources were replaced 
in the first self-calibration loop with point sources with the 
same centroid and flux density as the original source. Further 
self-calibration loops were done using the clean components 
in the traditional manner. Each day's observation was ini- 
tially self calibrated alone using 4 loops of self-calibration (2 
loops of phase calibration and 2 amplitude and phase cal- 
ibration loops). At this point the bright continuum source 
J104420-011146 was subtracted from the u-v data, and the 
data were flagged to exclude any visibility residuals that 
exceeded a threshold set by the system noise statistics. The 
bright continuum source was then added back into that day's 
data. 

After this detailed editing process, all the u-v data were 
combined and self-calibration was repeated on the entire 
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combined data set to ensure consistency between the dif- 
ferent days. A large continuum image of the entire field was 
then made. 

The AIPS imaging routine assumes that the sky is flat 
(it ignores the vertical 'w' term) which is appropriate only 
for small fields. In order to reduce the distortion over our 
large field it was necessary to break the image up into 
16 facets (see Figure 1). Each of these facets was 11.95' 
square (1024 pixels) and each overlapped by 1.2'. The final 
total field size was 44.20' on each side (i.e. large enough to 
include all the Fujita galaxies) . Clean boxes were put around 
all the significant continuum sources during the imaging pro- 
cess. In this initial continuum image, faint sidelobes from the 
brightest radio source J104420-011146 remained visible (see 
Figure 2). To reduce this, all the weaker continuum sources 
were subtracted from the u-v data, and the data were once 
again self-calibrated using clean components of the bright 
source alone. This new model of the bright source was sub- 
tracted from the original u-v data, and the weak sources 
(along with the appropriate calibration for these sources) 
restored. The final continuum image of the field without the 
bright source was then made. 

The final image has a resolution of ~2.9". The RMS at 
the centre of the final continuum image is 15 /J,3y. This is, to 
our knowledge, the most sensitive image achieved so far at 
the GMRT. The astrometry of the bright radio continuum 
sources in the data were checked against their positions in 
the VLA FIRST survey. A small offset of (-0.79 ± 0.30)" in 
right ascension and ( — 1.42±0.28)" in declination was found. 
This offset is believed to have been introduced to the data 
during the self-calibrating process though the exact origin of 
the offset is unclear. After correcting for this, several objects 
could be identified with aligned optical and radio continuum 
components. 

The continuum emission was then subtracted from the 
original spectral line data using the AIPS task UVSUB. 
While this removed most of the emission, small residuals of 
the brighter sources remained. To remove these final traces, 
a linear fit to the continuum across frequency was subtracted 
in the image plane using the AIPS task IMLIN. This intro- 
duces a small bias, since any line emission from the galaxies 
would be included in the fit. The correction for this bias is 
described in Section 4.1. In the final data cubes, the median 
RMS was 123 /iJy per channel in the lower sideband and 
136 /iJy per channel in the upper sideband. 



3 THE RADIO CONTINUUM AND Ha 
LUMINOSITY CORRELATION 

Synchrotron radio emission is generated in areas of active 
star formation from relativistic electrons accelerated in su- 
pernova remnants. Radio emission at 1.4 GHz generated by 
this process has been used as a measure of the star for- 
mation rate for galaxies, and it correlates well with other 
star formation indicators (Sullivan et al. 2001). In our ra- 
dio continuum image, it is possible to measure this emission 
for some of the Fujita galaxies. Since these galaxies are at 
a redshift of z = 0.24 and we observed at 1150 MHz, the 
rest frequency of their radio continuum emission in our data 
is 1.4 GHz. This means that we can directly compare the 
relationship in star forming galaxies between Ha line emis- 
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Figure 3. Ha emission line luminosity plotted against the 
restframe 1.4 GHz radio continuum luminosity for galaxies at 
z = 0.24. The line plotted through the data is the observed re- 
lationship found at z = (Figure 2 of Sullivan et al. 2001). The 
4 points in the 'individual galaxies' region are the brightest Ha- 
emitting galaxies. The triangle is a merging/interacting system 
of two galaxies. The square in the 'bright L(Ha)' region is the 
combined signal from galaxies with bright Ha emission. The two 
radio continuum 2cr upper limits in the 'medium L(Ha)' and 
'faint L(Ha)' regions are from the combined signal of subsam- 
ples of Ha-emitting galaxies. The star formation rate shown on 
the right axis is derived using the Kennicutt (1998) Ha to SFR 
relationship. 



sion luminosity and 1.4 GHz radio continuum luminosity at 
z = 0.24 and z ~ 0. Our measurements are shown in Fig- 
ure 3. Also shown in Figure 3 is the linear relationship found 
by Sullivan et al. (2001) from comparing the Ha emission 
against 1.4 GHz luminosity for 17 galaxies at z ~ 0. The 
uncertainties and the standard errors in the means for the 
Ha measurements in Figure 3 are all smaller than the size 
of the points. 

The integral flux densities of the radio continuum 
sources were measured using the AIPS task JMFIT which 
fits a 2D Gaussian to the source. The radio continuum for the 
4 brightest Ha-emitting galaxies each have individual mea- 
surable radio continuum fluxes (from ~100 to ~300 /iJy). 
These are the points shown in the 'individual galaxies' re- 
gion of Figure 3. The three galaxies with circular points show 
reasonable agreement with the Sullivan et al. line. The tri- 
angle point is a merging/interacting system of two galaxies 
within 3" (~11 kpc) of each other (see Figure 4). The galax- 
ies were observed to be at the same redshift of z = 0.2470 and 
have overlapping optical profiles. Both galaxies show strong 
Ha emission but only the galaxy with brighter Ha emission 
shows any signs of radio emission. Interestingly, this galaxy, 
although brighter in Ha, appears to be less massive than 
its companion (it is ~5 times more luminous in Ha than 
its companion but is ~0.6 times as luminous in the z' band 
filter). The reason why this galaxy's radio continuum lumi- 
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Figure 4. The top figure shows the Ho emission greyscale and ra- 
dio continuum contours of the pair of merging/interacting galax- 
ies mentioned in the text. The same pair of galaxies are shown in 
the bottom figure with the same radio continuum contours and 
the 'Iz continuum' (the 'Iz continuum' is the combination of the 
Suprime-Cam Iq and z' images, the continuum level around the 
Ha line; see Fujita et al. 2003 for details). The radio continuum 
contour levels are -40, -20, 0, 20, 40, 60, 80, 100 ^Jy- The images 
are both 7" on a side. 



nosity lies away from the Sullivan et al. line may have to do 
with it being a merging system. 

At lower Ha luminosities the radio continuum signal 
from multiple galaxies has been combined to improve the 
signal to noise ratio. Of the remaining 344 Fujita galaxies, 
4 have been excluded as they have contamination by bright 
radio continuum sources within 10" that do not appear to be 
connected to the galaxy. Another 2 galaxies were excluded 
because they lie within regions of artifacts in the synthesis 
images. Finally, we also exclude a galaxy with a very strong 
radio continuum flux (2336 ± 67 /iJy, which is brighter than 
any other Fujita galaxy by an order of magnitude) but has a 
relatively low Ha flux (39.82±0.10 log ergs" 1 ). The source 
of such a strong radio flux is more likely to be an active 



galactic nucleus than star formation. Unfortunately we do 
not have confirmation of the redshift for this galaxy. 

The galaxies have been divided up into three sub- 
samples for the purpose of combining the radio contin- 
uum. In these subsamples are: 45 bright galaxies with 
40.9 < L(Ha) 41.69 log ergs -1 , 55 medium galaxies 
with 40.4 < L(Ha) < 40.9 and 236 faint galaxies with 
39.65 < L(Ha) sC 40.4. Nearly all of the bright L(Ha) galax- 
ies and most of the medium L(Ha) galaxies have optical 
redshifts, but the vast majority of the faint L(Ha) galaxies 
do not. The radio continuum signal of the galaxies are com- 
bined using a weighted average, with the weight being the 
individual noise level for that galaxy. The further a galaxy 
is from the centre of the GMRT observations the greater the 
beam correction to its flux density; this results in a higher 
noise level and thus a lower weight in the sum. 

The square in Figure 3 is the average signal for the 
bright L(Ha) galaxies. This value shows good agreement 
with the Sullivan et al. (2001) line. For the medium L(Ha) 
and faint L(Ha) galaxies there was no detection of radio 
continuum flux, meaning that only upper limits can be de- 
termined. The 2cr upper limits for these regions are plotted 
on Figure 3. When determining these upper limits, it was 
assumed that the galaxies were not resolved by the GMRT 
synthesised beam. Both of these upper limits lie just to the 
left of the Sullivan et al. line, suggesting that the linear rela- 
tionship might not hold at low luminosities. Indeed, Sullivan 
et al. state that there is some evidence that the Ha-radio 
continuum relationship may not be perfectly linear for galax- 
ies with low star formation rates. They suggest that a frac- 
tion of the cosmic rays accelerated in supernovae remnants 
may escape from the lower mass galaxies. This would reduce 
the radio continuum emission relative to the galaxy's Ha lu- 
minosity. This effect would explain the positions of the two 
lower limits without requiring a change in the star formation 
mechanism at z = 0.24. 

Overall, Figure 3 shows that at z = 0.24 the relation- 
ship between radio continuum luminosity at 1.4 GHz and 
Ha emission luminosity for galaxies is consistent with that 
found at z = 0. This suggests that there has been no signif- 
icant change in the star formation mechanism over the last 
3 Gyrs. 



4 HI RESULTS 

4.1 The HI 21 cm Emission Signal 

To measure the coadded HI signal of the Fujita galaxies, we 
first use the optical redshifts to determine the expected fre- 
quency of the 21 cm HI emission. Of the 154 optical galaxy 
redshifts obtained, 121 lie within the usable regions of the 
radio data cube (see Figure 5). The reasons the remain- 
ing 33 galaxies are unusable are either: (i) their redshifted 
HI frequency lies outside the frequency range covered by 
the GMRT data, (ii) their redshift lies close to the bound- 
ary of the two sidebands where the data quality is poor, 
or (iii) their redshift lies within 5 channels of strong radio 
interference in the lower sideband (this strong radio inter- 
ference is at ~1137 MHz, covering ~0.625 MHz). For the 
purposes of defining a spectral window for coadding the HI 
signal, we allow for an average HI profile for the galaxies 
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Figure 5. The distribution of the measured redshifts for the Fu- 
jita galaxies. There are 154 redshifts in total. The dashed line is 
the shape of the Suprime-Cam narrow-band filter NB816 after 
converting the filter wavelength to the equivalent redshift of the 
Ha line. The vertical dotted lines are the GMRT frequency limits 
converted to HI redshift; the central line is the boundary between 
the upper and lower sidebands of the GMRT data. The histogram 
shaded with the unbroken line shows the distribution of the 121 
redshifts usable for HI coadding. The histogram shaded with the 
broken line shows the distribution of the 33 unusable redshifts 
(see text for details). 

of the order of ~300 kms" 1 (Doyle et al. 2005). A velocity 
window width of 500 kms" 1 should contain all the signif- 
icant HI emission after factoring in the uncertainty in the 
optical redshifts (~70 kms -1 ). 

In our data, the HI signals from individual galaxies are 
mostly below reasonable detection limits. Using the HI mass 
function at z ~ (Zwaan et al. 2005), we estimate that 
there is only a ~15% chance of a direct detection (5a) of 
the HI flux from an individual galaxy within the volume 
probed by our GMRT observations. There is probably more 
HI mass in galaxies at z = 0.24 than at z ~ 0, making a 
direct detection more likely, however detecting significant 
numbers of galaxies is still unlikely. For this reason we have 
coadded the HI signal from multiple galaxies to dramatically 
increase the signal to noise ratio of the measurement of the 
HI content of the sample. 

As galaxies are not individually detectable in our data, 
it is impossible to take into account any structure their HI 
gas may have when determining their HI flux density. How- 
ever, the peak specific intensity of an unresolved source in 
a radio synthesis image is equal to the flux density of that 
source. This means, that if a galaxy is unresolved by the 
GMRT beam, we can take the value of the specific intensity 
at the optical position of the galaxy as a measure of the 
galaxy's total HI flux density. 

The GMRT synthesised beam has a FWHM ~2.9" 
which corresponds to ~11 kpc (1" corresponds to 3.8 kpc at 
z = 0.24). All but the smallest of the galaxies in our sample 
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Figure 6. This figure shows the histogram of the estimated 
effective HI diameter of the Fujita galaxies. This was derived 
from the diameter of the galaxy at the optical B magnitude 
25 mag arcsec -2 isophote and the correlation from Broeils & Rhee 
(1997). The histogram shaded with the dashed line shows the dis- 
tribution for all 348 Fujita galaxies; the histogram shaded with 
the solid line shows the distribution for galaxies with measured 
optical redshifts. The vertical dashed lines show the boundaries 
used in sorting galaxies into the three different synthesised beam 
subsamples with beam FWHM of 2.9" (11 kpc), 5.3" (15 kpc) 
and 8.0" (20 kpc). 

are at least partially resolved by this beam size. To solve this 
problem the radio data were smoothed to larger beam sizes 
of 5.3" (~20 kpc) and 8.0" (~30 kpc). The larger galaxies 
should be unresolved in these larger beam size data. Gaus- 
sian smoothing in the image plane is equivalent to tapering 
(multiplying by a Gaussian) in the u-v plane. This effectively 
reduces the weight of the outer antennas in the image plane, 
and hence decreases the signal to noise ratio for unresolved 
sources. However the measured HI flux increases for galaxies 
that are now unresolved in the smoothed data. 

A rough estimate of the HI size of each individual galaxy 
was made using a correlation between optical and HI size 
found by Broeils & Rhee (1997). From their data, we derive 
a relationship between the diameter of the optical isophote 
corresponding to a surface brightness of 25 mag arcsec -2 
in B (D25) and the diameter within which half the HI mass 
of the galaxy is contained (D e g). This relationship is: 

log(Deff) = (0.978 ± 0.035) log(D^) + (0.041 ± 0.046) (1) 

where D^g 1 and D e d are measured in kpc. 

The IRAF task ellipse was used to measure the 
25 mag arcsec -2 isophotal level in the Subaru B images, 
taking into account the cosmological surface brightness dim- 
ming and a k-correction taken from Norberg et al. (2002). 
From this an estimate of the HI diameter is made. 

If a galaxy has 2 D c ff < D m id, where D m id is the mid- 
point between two beam smoothing levels, then it was as- 
sumed to be unresolved in the smaller beam data. Figure 6 
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Figure 7. The average HI galaxy spectrum created from coadding 
the signal of all 121 galaxies with known optical redshifts. The 
top spectrum has no smoothing or binning and has a velocity 
step size of 32.6 kms" 1 . The bottom spectrum has been binned 
to ~500 kms -1 . This is the velocity width that the combined HI 
signal of the galaxies is expected to span. For both spectra the 
1<t error is shown as dashed lines above and below zero. 



shows the distribution of the estimated effective HI diame- 
ter and the grouping of the galaxies into the three different 
beam sized radio data sets. 

When making the final data cube, a linear fit to the 
spectrum through each sky pixel was subtracted to remove 
the residues left from the continuum sources (see discussion 
on IMLIN previously). Any spectral line features are in- 
cluded in the calculation of this linear fit across frequency. 
Such features would create an over-estimation of the con- 
tinuum, creating a bias in the measured HI spectrum. To 
correct for this effect, each galaxy HI spectrum has a new 
linear fit made across all channels except those channels cor- 
responding to a 500 kms -1 velocity width around the galaxy 
at its redshift (these should be the channels that contain any 
HI signal). This new fit is then subtracted from the data, 
correcting for any bias created by IMLIN. This correction 
to the HI flux for the combined galaxy spectrum increases 
the measured flux by ~25%. 

When coadding the separate HI spectra, each measured 
flux value carries a statistical weight that is inversely pro- 
portional to the square of its uncertainty (the RMS noise 
level). This noise value is calculated from the known RMS 
per frequency channel in the radio data cube, factoring in the 
increase in the noise from the beam correction for galaxies 
away from the centre of the GMRT beam. As we are probing 
far from the GMRT beam centre, uncertainties in the beam 
shape could be a significant factor. However as the coadded 
HI signal from individual galaxies is weighted by their noise, 
any errors in the assumed Gaussian shape in the outskirts 
of the beam are heavily attenuated by the weight factor. 

The weighted average HI spectrum from coadding all 
121 galaxies can be seen in Figure 7. To measure the error 



in the HI spectrum, a series of artificial galaxies with ran- 
dom positions and HI redshifts were used to create coadded 
spectra. From many such artificial spectra a good estimate of 
the noise level in our final real spectrum can be determined. 
As seen in Figure 7, the noise level increases with increasing 
velocity offset from the centre of the coadded spectrum. This 
is because some galaxies lie at redshifts near the edges of the 
data cube and when adding the spectra of these galaxies to 
the total, there is no data for channels that lie off the edge 
of the data cube. These channels with no data are given zero 
weight in the coadded sum resulting in a higher noise level 
at these velocities in the final coadded spectrum. 

In Figure 7 the total HI flux density within a central 
velocity width of 500 kms -1 is 8.9 ± 3.4 mjykms -1 , a sig- 
nal to noise ratio of 2.6. The level of this flux is not highly 
dependent on the choice of the velocity width; choosing var- 
ious velocity widths from ~400 to ~600 kms -1 gives total 
fluxes that are similar within the errors. 

HI emission flux can be converted to the mass of atomic 
hydrogen that produced the signal by the following relation: 



Mi 



236 



S v 

( 1 + z ) V mjy 



Mpc 



AV 

km s~ 



(2) 



where S v is the HI emission flux averaged across the velocity 
width AV and dL is the luminosity distance to the source. 
This equation assumes that the cloud of atomic hydrogen gas 
has a spin temperature well above the cosmic background 
temperature, that collisional excitation is the dominant pro- 
cess, and that the cloud is optically thin (Wieringa et al. 
1992). No correction for HI self absorption has been made. 
HI self-absorption may cause an underestimation of the HI 
flux by as much as 15%. However this value is extremely un- 
certain (Zwaan et al. 1997). For our HI mass calculation we 
used AV = 500 kms -1 , mean z = 0.235 and d L = 1177 Mpc. 
This gives an average HI galaxy mass for all 121 galaxies of 
(2.26 ±0.90) x 10 9 M Q (Note: Mhi = 6.3 x 10 9 M at z = 0, 
Zwaan et al. 2005). 



4.2 The Cosmic Neutral Gas Density, f2 gas 

The conversion of the HI mass of the galaxies to the cosmic 
density of neutral gas requires a measure of the number den- 
sity of the galaxies. This information is available from the 
Fujita et al. (2003) Ha luminosity function at z = 0.24. 

The narrow-band filter used to select the Fujita galax- 
ies is not square but has an almost Gaussian-like profile (see 
Figure 5) . This means that galaxies with strong Ha emission 
can be detected in the wings of the filter but will have their 
measured Ha luminosity underestimated. There is a strong 
luminosity bias in the spectroscopic redshifts we obtained 
due to this effect. Bright galaxies that have Ha emission 
lines in the wings of the filter will be easier to obtain red- 
shifts for than faint galaxies at the centre of the filter. The 
distribution of redshifts appears to be slightly bias towards 
lower redshifts. This small bias is probably caused by the 
[Nil] 6584 emission line contributing more to the narrow- 
band flux at the lower redshifts, giving rise to more detec- 
tions. Near z ~ 0.225, the emission from the [NII]6584 line 
and the [SII] doublet both come through simultaneously at 
opposite ends of the narrow-band filter. This could explain 
the tail at low redshift seen in the redshift distribution. 
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Galaxy 
Subsample 


Number 

of 
galaxies 


Limits 

(log erg s _1 ) 


Average 
L(Ha) 

(log erg s _1 ) 


Average 
SFR 

(M yr- 1 ) 


Average 
HI Mass 

(10 9 M ) 


Average 
Volume per 
galaxy (Mpc 3 ) 


Subsample 

figas 

(io- 3 ) 


all galaxies 


121 


— 


41.0 


0.79 


2.26 ± 0.90 


— 


— 


bright L(Ha) 


42 


L(Ha) > 40.9 


41.4 


1.98 


4.1 ± 2.2 


105 ± 38 


0.38 ± 0.24 


medium L(Ha) 


42 


40.4 < L(Ha) ^ 40.9 


40.6 


0.31 


2.9 ± 1.5 


76 ± 28 


0.37 ±0.23 


faint L(Ha) 


37 


40.0 < L(Ha) < 40.4 


40.2 


0.13 


0.9 ± 1.4 


53 ±27 


0.16 ± 0.26 



Table 1. The table shows the properties of all the galaxies coadded together and the three Ha luminosity subsamples. For each subsample 
the average Ha luminosity, average star formation rate (SFR) and average HI mass is listed. Also listed is the average volume per galaxy 
which is the average volume one would need to probe in order find a single galaxy belonging to that luminosity subsample. The reciprocal 
of this volume is the number density of such galaxies per Mpc 3 . The cosmological neutral gas density for each galaxy subsample is listed 
in the last column. Combining these values gives the total cosmic neutral gas density. The star formation rate listed is derived using the 
Kennicutt (1998) Ha to SFR relationship. 



Fujita et al. (2003) applied a statistical correction for 
the effect of the filter shape when creating their Ha luminos- 
ity function. This could be done because they have reason- 
ably well-understood selection criteria for Ha emission de- 
tection. This is not the case for the Fujita galaxies for which 
we obtained redshifts. Different galaxies were observed for 
different amounts of time, as it was not possible to configure 
the fibres in 2dF on all the galaxy targets simultaneously. 
Some galaxies also had small errors in their astrometry that 
reduced the amount of light down the optic fibre obtained 
for that galaxy (these astrometric errors have since been 
corrected). This uncertainty in the detection success rate 
for the galaxy redshifts is why we are using the Fujita et al. 
Ha luminosity function to determine the number density of 
the galaxies. 

For calculating the number density of the galaxies, the 
Fujita galaxies with redshifts have been divided into three 
subsamples based on their Ha luminosity. In order to di- 
vide the galaxies accurately into subsamples, it is neces- 
sary to correct the Ha luminosity for the effect of the fil- 
ter shape. This is done in order to separate genuine faint 
galaxies from brighter galaxies seen through the wings of 
the narrow-band filter. Without this correction, the faintest 
Ha emission galaxies are measured to have more atomic hy- 
drogen gas than galaxies with stronger Ha emission. 

An estimate of the necessary Ha luminosity correction 
can be made using the precise redshifts of the galaxies and 
the known shape of the narrow-band filter with wavelength 3 . 
The measured Ha luminosity is adjusted by the relative de- 
crease in the filter transmission at the redshifted wavelength 
of the Ha line compared to the filter centre 4 . 



3 National Astronomical Observatory of Japan website 
http: / / subarutclcscope.org/ Obscrving/Instruments/SCam / 

4 This Ha correction was not used in the comparison of the Ha 
and radio continuum observations in Section 3. Corrections can- 
not be done for all the galaxies included in that analysis as not 
all have redshifts. The correction is only necessary in the work 
above because of the bias in the galaxies for which we obtained 
spectroscopic redshifts. The much larger sample of faint galaxies 
in Section 3, most of which do not have redshifts, is dominated by 
genuine faint galaxies rather than bright galaxies shining through 
the edge of the narrow-band filter. It is these bright galaxies that 
need the correction. The number of these bright galaxies mis- 



Using the corrected Ha luminosities the galaxies are 
separated into three subsamples: bright, medium and faint. 
The details for each of the Ha subsamples can be seen in 
Table 1. The average HI mass for the galaxies in each sub- 
sample is measured as described in Section 4.1. 

An estimate of the average number density for galaxies 
belonging to each subsample is determined using the Fujita 
et al. (2003) luminosity function. This is done by integrat- 
ing the Ha luminosity function over the Ha limits for each 
subsample to obtain the number density of such galaxies. 
The galaxy number density for each subsample is expressed 
in Table 1 as the average volume one would need to probe 
in order to find a single galaxy belonging to that Ha lu- 
minosity subsample (the reciprocal of the number density) . 
The errors in the calculated volumes are derived by consid- 
ering the effect of the errors in the Fujita et al. Schechter 
luminosity function parameters. 

The cosmic density of HI gas for each subsample can be 
calculated using: 



where p hi is the cosmic density of HI gas, Mm is the aver- 
age HI mass per galaxy and V ga i is the average volume per 
galaxy. To compare with other values in the literature, it is 
necessary to convert this value of the HI gas density to the 
density of neutral gas. This is done by adding a correction 
for the neutral helium content, which is assumed to be 24% 
by mass of the total neutral gas. We convert this neutral gas 
density to units of the total cosmic mass density using: 

_ 8 TV G p gas 

"gas — gjj2 W 

After combining the neutral gas density values for each sub- 
sample we obtain a total cosmic density of neutral gas in star 
forming galaxies at z = 0.24 of f2 gas = (0.91 ± 0.42) x 10~ 3 . 
This value is shown in comparison to other literature values 
in Figure 8. 

Figure 8 shows the neutral gas density of the universe 

classified as faint galaxies should be small compared to the large 
number of actual faint galaxies in that analysis, so that any cor- 
rection should not be statistically significant. 
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Figure 8. The neutral gas density of the universe as a function of redshift (on the left) and look back time (on the right). The small 
triangle at z = is the HIPASS 21 cm emission measurement from Zwaan et al. (2005). The filled circles are damped Lya measurements 
from Prochaska et al. (2005). The open circles are damped Lya measurements from Rao et al. (2006) using HST. The large triangle 
at z = 0.24 is our HI 21 cm emission measurement made using the GMRT. All results have been corrected to the same cosmology and 
include an adjustment for the neutral helium content. 



as a function of both redshift and look-backtime. There is 
a large uncertainty in the cosmic neutral gas density in the 
redshift range z = 0.1-1.5, which corresponds to two thirds 
of the age of the universe. Our measured value of £2 gas is com- 
parable to that from previous damped Lya measurements 
at intermediate redshifts (Rao et al. 2006). While the trend 
of the measurements is to indicate a sharp rise in HI content 
with redshift, both studies suffer from large statistical un- 
certainty. However, our new measurement has the advantage 
that it applies to a narrow range of redshifts z ~ 0.24. 

Strictly speaking, our measured cosmic neutral gas den- 
sity is only a lower limit, as it includes only the gas from 
galaxies that show star formation. The star formation-HI 
mass correlation found in low-redshift Hl-selected galaxies 
(Doyle & Drinkwater 2006), shows that galaxies with little 
or no star formation contain no significant quantities of neu- 
tral gas. This relationship appears to hold at z = 0.24 (see 
Section 4.3). Since our sample includes galaxies with Ha lu- 
minosities at ~1% of L* no significant reservoirs of HI gas 
should be missed. 

In this work no attempt has been made to separate the 
Ha emitting galaxies into those that are genuine active star 
forming galaxies from those that are actually active galactic 
nuclei (AGN). When calculating their star- formation den- 
sity for this field at z = 0.24, Fujita et al. (2003) included 
a 15% correction for the luminosity contribution of AGN. 
(This correction came from Pascual et al. 2001). No simi- 
lar correction exists for the atomic hydrogen gas content of 
AGN. However, the ratio of HI mass to Ha luminosity for 
AGN is likely to be lower than for star forming galaxies, as 
AGN preferentially reside in early-type galaxies (Schade et 
al. 2001). Including AGN in our sample may cause a small 
underestimation in the density of neutral gas for star form- 



ing galaxies, however it will bring our result closer to the 
true total cosmic neutral gas density of all galaxies by in- 
cluding another galaxy type in the measurement. 

The volume probed by Fujita et al. (2003) for their Ha 
luminoisty function (and used here for the calculation of the 
number density of galaxies) is ~4000 Mpc and is a rectan- 
gular strip with a depth ~110 Mpc. The cosmic variance in 
the number of galaxies in a field of this volume and shape 
at z = 0.24 is approximately 40%. This estimate is based on 
a calculation using the Eisenstein & Hu (1998) fitting func- 
tion for the matter power spectrum functon and is confirmed 
by an examination of the galaxy distribution in the Millen- 
nium Simulation at z = 0.24 (Springel et al. 2005; Croton 
et al. 2006). Treating this as a random error (rather than a 
systematic offset), the final error in f2 gas would increase by 
~15%. However, the star formation rate measured for this 
field places it close to the mean, or perhaps slightly higher, 
compared to other values in the literature (see Section 2.1). 
This suggests that this field has close to the average num- 
ber of galaxies despite the large possible cosmic variation, 
and that the measured neutral gas density measured in this 
field is close to, or slightly higher than, the true cosmic gas 
density, f2 gas . 



4.3 The Star Formation Rate-HI Mass Galaxy 
Correlation 

In the local universe there is a known correlation between 
the star formation rate in a galaxy and the mass of HI in that 
galaxy. For a recent study using Hl-selected galaxies from 
HIPASS, see Doyle & Drinkwater (2006). In Table 1 we have 
listed the average HI galaxy mass for the different Ha lu- 
minosity subsamples. For each subsample, one can derive 
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Figure 9. The average galaxy atomic hydrogen gas mass plotted 
against the average galaxy star formation rate. The circular points 
are the average values for the bright, medium and faint L(H«) 
subsamples (see Table 1). The triangle is the value for the average 
of all galaxies. The line is from Doyle & Drinkwater (2006) for 
z ~ galaxies. 



5 CONCLUSION 

We have demonstrated that is possible to measure, in rea- 
sonable observing times, the average amount of atomic hy- 
drogen gas in galaxies from their HI 21 cm emission at in- 
termediate redshifts. This is done by coadding the HI signal 
from multiple galaxies with known positions and optical red- 
shifts. Our measurement of HI 21 cm emission from galaxies 
corresponds to a cosmic density of neutral gas that is consis- 
tent with previous measurements of the neutral gas density 
in damped Lya measurements at a similar redshift. 

The relationship between Ha luminosity and the rest- 
frame 1.4 GHz radio continuum emission in star-forming 
galaxies at z = 0.24 has been shown to be consistent with 
the correlation found at z = 0, as is the the relationship 
between galaxy star formation rate and galaxy HI mass. 
These two results suggest that the process of star formation 
in field galaxies is not significantly different 3 Gyr ago from 
the present day. 

However, the significantly higher global star formation 
rate density found at still higher redshifts does not seem to 
continue this trend. At z > 1, the damped Lya measure- 
ments of the neutral gas density and the star formation rate 
density differ by a factor of two from what one would expect 
from the star formation-HI mass correlation found at z = 0. 
This implies an enhanced efficiency of star formation for a 
given neutral gas density at high redshifts. 



the average star formation rate per galaxy from the average 
Ha luminosity and the SFR conversion of Kennicutt (1998). 
Our comparison of star formation rate and galaxy HI mass 
is shown in Figure 9. Plotted on this figure is the linear 
relationship from Figure 3 of Doyle & Drinkwater (2006), 
where they compared individual galaxies' HI masses from 
HIPASS to their star formation rate derived from IRAS in- 
frared data. The linear fit, taken directly from their figure, 
is: 
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log(HI Mass) = 0.59 log(SFR) + 9.55 (5) 

From Figure 9, it is clear that the values found at z = 0.24 
are consistent with the correlation found at z = 0. This 
suggests that there has been no significant change in the 
star formation mechanisms in field galaxies 3 Gyrs ago; that 
the same amount of atomic hydrogen gas in a galaxy gives 
the same measurable star formation rate. 

At z > 1 there is a factor ~10 increase in the star for- 
mation density compare to z = (Hopkins 2004). If the 
star formation-HI mass correlation as stated in equation 5 
were to hold true at these high redshifts, a factor of ~4 in- 
crease in neutral gas density would be required. However, the 
damped Lya measurements of Prochaska et al. (2005) only 
indicate that the neutral gas density rises by factor of ~2. 
This implies that star formation rate-HI mass relationship 
must have been different at these higher redshifts, indicating 
a change in the mechanism of star formation. The mecha- 
nism for converting gas to stars must have been a factor of 
~2 more efficient at high redshifts. 
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